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We study the influence of particle shape on growth processes at the edges of evaporating drops. 
Aqueous suspensions of colloidal particles evaporate on glass slides, and convective flows during 
evaporation carry particles from drop center to drop edge, where they accumulate. The resulting 
particle deposits grow inhomogeneously from the edge in two-dimensions, and the deposition front, 
or growth line, varies spatio-temporally. Measurements of the fluctuations of the deposition front 
during evaporation enable us to identify distinct growth processes that depend strongly on parti- 
cle shape. Sphere deposition exhibits a classic Poisson like growth process; deposition of slightly 
anisotropic particles, however, belongs to the Kardar-Parisi-Zhang (KPZ) universality class, and de- 
position of highly anisotropic ellipsoids appears to belong to a third universality class, characterized 
by KPZ fluctuations in the presence of quenched disorder. 

PACS numbers: 61.43.Fs,64.70.kj,64.70.pv,82.70.Dd 
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Examples of surface and interfacial growth phenomena 
are diverse, ranging from the production of uniform coat- 
ings by vapor deposition of atoms onto a substrate [l|, 
to burning paper wherein the combustion front roughens 
as it spreads 0, to bacterial colonies whose bound- 
aries expand and fluctuate as bacteria replicate Q. The 
morphology of the resulting interfaces is a property that 
affects the macroscopic responses of such systems, and 
it is therefore desirable to relate interface morphology to 
the microscopic rules that govern growth 0, @] ■ To this 
end, simulations have directly compared a broad range of 
growth processes 0, @ and have found, for example, that 
the random deposition of repulsive particles is a Poisson 
process, while the random deposition of "sticky" parti- 
cles belongs to a different universality class that leads to 
different interface morphology. 



Besides discrete models, theoretical investigation of 
this problem has centered around continuum growth 
equations (e.g., Q). One interesting approach that uni- 
fied a large set of discrete simulations is based on the 
so-called Kardar-Parisi-Zhang (KPZ) equation d, H-12|. 
This nonlinear equation relates stochastic growth and in- 
terfacial growth fronts/lines/surfaces to diffusion and lo- 
cal lateral correlations; its solutions are known and be- 
long to the KPZ universality class (l3Wl(i|. The KPZ class 
presents a rare opportunity for connecting exact theoreti- 
cal predictions about nonequilibrium growth phenomena 
with experiment. However, to date only a few members 
of the KPZ class have been experimentally identified 0- 
This paucity of KPZ examples is due, in part, to 
the presence of quenched disorder and long-range interac- 
tions in experiment, as well as to limited statistics, which 
make growth process differences difficult to discern. In 
fact, broadly speaking, experimental confrontation of the 
microscopic rules explored by theory and simulation has 



been difficult. 

In this contribution we demonstrate that the rich 
nonequilibrium physics of evaporating colloidal drops 
provides an attractive experimental system for study of 
such growth processes and for testing theory and simula- 
tion predictions. Specifically, we investigate the growth 
of particle deposits from the edges of evaporating col- 
loidal drops. Particle deposition is observed by video 
microscopy at the single particle level. Aqueous suspen- 
sions of colloidal particles are allowed to evaporate on 
glass slides at constant temperature and humidity, and 
radial convective flows during evaporation carry particles 
from drop center to drop edge, where they accumulate 
(Fig. la) [18(. The resulting deposits of particles grow 



from the edge in two-dimensions, defining a deposition 
front, or growth line, that varies in space and time. Inter- 
estingly, these interfacial growth processes are strongly 
dependent on colloidal particle shape [lj|. Three dis- 
tinct growth processes were discovered in the evaporating 
colloidal suspensions by tuning particle shape-dependent 
capillary interactions and thus varying the microscopic 
rules of deposition. The substantial shape fluctuations 
of the growth line of spheres are readily explained via a 
Poisson like deposition process; slightly anisotropic par- 
ticles exhibit weaker fluctuations characteristic of KPZ 
class behavior, and very anisotropic ellipsoids exhibit be- 
havior consistent with the KPZ class in the presence of 
quenched disorder [20l - [22| . 

Our experiments employ water drops containing a 
suspension of micron-sized polystyrene spheres (Invitro- 
gen) stretched asymmetrically to different aspect ratios 
USE Hi- The spheres are 1.3 pm in diameter; all el- 
lipsoids are stretched from these same 1.3 pm spheres. 
We evaporate the drops on glass slides (Fisher Scien- 
tific) and study suspensions containing particles of the 
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FIG. 1: a. Cartoon depicting the deposition mechanism. A 
radially outward flow carries particles from the drop center 
to the drop edge, where they are deposited on the air-water 
interface, b-d. Binarized experimental images of deposits of 
spheres (e = 1.0) (b), slightly stretched particles (e = 1.2) 
(c), and ellipsoids (e = 2.5) (d), along with images of single 
particles. A label in (b) demonstrates our definition of height 
(h), i.e., distance from drop edge. The drop edges and the 
direction of the coffee-ring driving flow are indicated. The 
probed window size L is defined for L = 100 pm (b) and 
L = 25 pm (d). These scale bars also hold for (c). e-g. The 
deposit height profile (growth line), h, plotted as a function 
lateral position, x, at four different times, for e = 1.0, 1.2 and 
2.5 (e-g, respectively). 



same composition, but with different major-minor diam- 
eter aspect ratio (e), including spheres (e = 1.0), slightly 
anisotropic particles (e = 1.05, 1.1, 1.2), and ellipsoids 
(e = 1.5,2.5,3.5). We study volume fractions (4>) that 
vary from <j> = 10 -4 to 0.02. 

The experiments are reproducible across many 
droplets, enabling the accumulation of sufficient statis- 
tics to test continuum equation predictions of surface 
roughness scaling, and more. Further, the strong shape- 
based capillary attractions between particles on the air- 
water interface f25T - l3lj ] permit us to establish relation- 
ships between particle interaction and interfacial growth 
processes. 

Qualitative differences between deposits crafted from 
particles with different aspect ratio, e, are readily appar- 
ent (Fig. lb-d). Deposits of spheres are densely packed 
(Fig. lb); deposits of slightly anisotropic particles are 
loosely packed (Fig. lc); deposits of very anisotropic el- 
lipsoids form open networks and "empty" structures (Fig. 
Id). Deposits are characterized by their height, h, mea- 
sured from the three-phase contact line (i.e., their pen- 
etration from the drop edge, measured radially inward) 
(see Fig. lb). The deposit morphology is well described 
by the contour of its height profile (growth line). Almost 
all of the deposits described in this letter are monolayers 
of particles adsorbed on the air- water interface. The sole 
exceptions are the very largest sphere deposits (e = 1.0) 
discussed below. 

Deposit height varies spatially and temporally 
(h(x,t)). For example, the mean height, h, increases in 
time (see [32]). Spatial variation of h is quantified by 
the standard deviation of h, a quantity referred to as the 



FIG. 2: Deposit width (w), i.e., the standard deviation 
of deposit height (h), plotted versus h for spheres (e = 1.0, 
squares), slightly anisotropic particles (e = 1.2, circles), and 
ellipsoids (e = 2.5, 3.5, triangles). Different colored dots indi- 
cate data derived from different experiments. The data col- 
lapse onto three unique trend lines based on e. 



width, w. The growth of local fluctuations in h produce 
an increase in w over time (32j . Many growth processes 
exhibit a self-affine structure well-described by Family- 
Vicsek scaling. For example, w oc L a for small L, where 
L is the horizontal (lateral) size of the window in which 
w is calculated. Similarly, the width can exhibit power- 
law growth over time, i.e., w « r. Here a is called the 
roughness exponent and /3 is called the growth exponent 
For simple Poisson processes, /3 — 1/2, but w does 
not depend on L, so a is poorly defined. For members of 
the KPZ universality class, (3 — 1/3 and a = 1/2. 

In exploring growth exponents, the use of real time, 
t, presented us with technical problems. During evap- 
oration, for example, the radially outward flow rate in- 
creases over time [181 ] , preventing a controlled power law 
measurement. Further, the start time (t = 0) is ambigu- 
ously defined; it could be when the drop is placed on the 
substrate, when the drop stops spreading, or when the 
first particle is deposited. Finally, the time-range during 
which deposits can form, i.e., before strong surface flows 
or particle aggregates artificially modify w, is often not 
long enough to extract meaningful power laws. By con- 
trast, the use of average deposite height, h, in place of 
time t averts many of these technical issues, h is linear 
with t in all experiments, with a coefficient of determi- 
nation R 2 > 0.9 in all cases. Importantly, use of h skirts 
the issue of increasing/decreasing flow rates, resolves am- 
biguities in defining t — (h = is well defined), and 
enables us to perform many measurements with differ- 
ent drops, or on different parts of the same drop, all of 
which can be combined (Fig. 2c). Thus, we are able to 
obtain sufficient statistics to provide data over two orders 
of magnitude in h. 

Three qualitatively and quantitatively different growth 
rate regimes are readily identifiable in the results (Fig. 
2). For spheres (e = 1.0) to oc h 0AS ^; for slightly 
anisotropic particles (e = 1.2) w oc ft°- 37 ( 4 ); for ellip- 
soids with e = 2.5 and e — 3.5, the data collapse onto 
a single curve with w oc /i - 68 ( 4 ) (the final digit uncer- 
tainty represents the standard error of the fit combined 
with measurement uncertainty). Three distinct growth 
processes are readily apparent. Spheres pack densely at 
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the drop edge, but their growth deposit appears spatially 
uncorrelated; thus we observe a significant increase in w 
over two decades in h. Slightly anisotropic particles start 
with a relatively large w, which then increases slowly with 
h. For the anisotropic ellipsoids, w increases rapidly with 
ft, as dense regions grow at the expense of sparse regions, 
which remain sparse. 

The measured growth exponent for slightly anisotropic 
particles (e = 1.2), /3 = 0.37(4), is consistent with the 
KPZ universality class, which predicts ,9 = 1/3 and 
a = 1/2. We measure a based on the finite size scal- 
ing of w, as w depends on the observation lengthscale, 
L. For small values of L, w oc L a . (Note, at larger L, 
w saturates or crosses over to a weaker power-law de- 
pendence.) The power-law dependence is measured over 
two decades in L within the small L regime (Fig. 3a). 
The best fit yields a = 0.51(5). Both scaling exponents, 
P = 0.37(4) and a — 0.51(5), are within experimental 
uncertainty of predicted values for the KPZ universality 
class. 
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FIG. 3: a. Width, w, plotted versus probed lengthscale, 
L, for a drop containing slightly anisotropic particles e = 1.2 
with height h = 13 jim. The dashed line represents the best 
power law fit for L < 15 pm. The observed power law scaling 
is consistent with the KPZ universality class, b. Skewness 
(71) and kurtosis (72) of the h distribution for particles with 
anisotropy e = 1.2 is plotted versus h. Solid lines represent 
values of 71 and 72 for Gaussian unitary ensemble matrices. 
Thus, the skewness and kurtosis imply that the distribution 
of height is consistent with the KPZ universality class. 



Next, we characterize the distribution of ft, 
in samples containing particles with anisotropy 
e = 1.2. In particular, the skewness, 71 (ft) = 
(1/L J x=Q L ((h(x,h) — h) 3 / 2 )/w(h) 3 )dx, and excess kur- 
tosis 7 2(ft) = (1/L J x=0 L ((h(x,h) - h) A )/w(h) 2 )dx - 3, 
are calculated; here h(x, ft) is the height at position x, 
at a time with mean height ft. For members of the KPZ 
class, the distribution of ft depends on the shape of the 
interface ll|. For curved interfaces (e.g., the circular 
three-phase contact line of a sessile drop), a Gaussian 
unitary ensemble (GUE) random matrix distribution 
is expected (71 = 0.22 and 72 = 0.09), while for flat 



interfaces, a Gaussian orthogonal ensemble distribution 
is expected (71 = 0.29 and 72 = 0.16) [n[. In our 
experiments, the height distribution levels out at ft w 15 
um where 71 = 0.20(6) and 72 = 0.10(3) (Fig. 3b). Thus, 



the distribution of height is consistent with predictions 
for a KPZ process with curved interface. 

Finally, we note that particles with anisotropy e = 1.2 
adsorb on the air-water interface and slightly deform the 
air- water interface [25| . These deformations, in turn, in- 
duce a relatively weak capillary interparticle attraction. 
The attraction is significant at short-range, enabling par- 
ticles to "stick" to each other once they reach the air- 
water interface. This behavior is similar to simulations 
of so-called "ballistic deposition" [(| which gives rise to 
KPZ dynamics. 

For spheres, the measured growth exponent ft = 
0.48(4) (for ft < 10 pm) is consistent with a Poisson pro- 
cess (raining particles) which might be expected, since 
spheres do not significantly deform the air-water inter- 
face, and thus do not induce the strong quadrupolar at- 
traction observed for ellipsoids. In this case, the distri- 
bution of heights is well fit by the Poisson distribution, 
N(h) — X exp(— X)/hl, where N(h) is the number of oc- 
currences of ft, and A is the mean of the distribution. 
The best fit A increases linearly with h {R 2 = 0.99) [12]. 
The standard deviation of a Poisson distribution is vft, 



and 71 = ft 



-1/2 



i.e., the skewness of the height distri- 



bution decreases as the width increases. In our experi- 
ments, 71 does not approach its asymptotic values until 
ft > 1.0 pm, so a power law fit could only be measured 
over one decade. Instead, to determine whether these 
data are consistent with Poisson distributions, we plot 
71 • w (Fig. 4a). Since w oc ft 05 , we expect 71 • w to be 
a constant when 71 oc w^ 1 m ft -05 . Experimentally, for 
e = 1.0, 71 oc ft, -0 5 , consistent with a Poisson process, 
and for e > 1.0, 71 • w is not consistent with a Poisson 
process. 

We next characterized the lateral correlations in w(h) 
for spheres. Unlike slightly anisotropic particle s, sp heres 
do not exhibit a distinct roughness exponent a 32( . The 
apparent absence of a characteristic roughness implies 
there is little spatial correlation in ft, i.e., regions with 
large ft may directly neighbor regions with small ftj which 
is again consistent with the Poisson distribution [32J . We 
note, however, that when the deposit is large (i.e., ft > 10 
pm), the deposit becomes multilayered and in three di- 
mensions, more contacts are required to stabilize a parti- 
cle, than in two dimensions. For h > 10 pm, the spheres 
are therefore able to find a local minimum in the height 
profile, rather than attaching to the first pair of particles 
they touch. The growth mechanism ceases to be a Pois- 
son like process of random deposition at this point, and 
it is better described as random deposition with surface 
diffusion [fjj]. 

For the very anisotropic ellipsoids, the measured 
growth exponent, j3 — 0.68(5), is consistent with a KPZ 
process in the presence of quenched disorder (KPZQ) 
Previous work found that if quenched dis- 



order prevents interfacial growth in a particular region, 
then a new universality class (the KPZQ class) was pro- 
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duced with /3 = 0.68 and a = 0.63 [2(J. The measured 
roughness exponent for our experiments is a = 0.61(2) 
(Fig. 4b). Thus the measured values of /3 = 0.68(5) and 
a = 0.61(2) are consistent with the KPZQ class. 

The latter observation was somewhat unexpected for 
us, because most KPZQ models are characterized by 
regions where growth is prevented, mixed with regions 
where growth is allowed ^HHll. Superficially, our ex- 
periments appear to have different conditions than those 
needed for KPZQ. However, the highly anisotropic ellip- 
soids induce strong capillary attraction on the air-water 
interface [25T - [3lj . which causes regions with many parti- 
cles to strongly attract additional particles. Even par- 
ticles that adsorb on the air-water interface in regions 
nearly void of particles are strongly attracted to particle 
rich regions, and eventually are deposited in particle rich 
regions. Thus, ellipsoids exhibit a colloidal "Matthew 



Effect" [34(. Surprisingly, this process is quite similar to 
the KPZQ scenario. The strong long-ranged capillary at- 
traction enhances growth in particle rich regions, which 
effectively prevents (or at least slows) growth in other 
regions (Fig. 4c). 

Finally, we explored a range of aspect ratios which fall 
between these exemplary cases. The dynamic scaling ex- 
ponent P is shown in Fig. 4d for many different values 
of e. Three regimes are readily identifiable. The dy- 
namic scaling exponent /3 starts at ~ 0.5 for spheres. As 
e increases, /3 abruptly decreases to ~ 0.35. When e is 
increased further, /3 increases to ~ 0.68. 

While quantitative differences between Poisson, KPZ, 
and KPZQ processes are provided by a and /3, qualitative 
differences can be summarized by performing a series of 
simple simulations 32] . Consider a system wherein par- 
ticles "rain down" vertically onto a ID surface divided 
into a series of columns (imagine Tetris with individual 
blocks) . A Poisson process can be modeled by randomly 
adding new particles without spatial or temporal correla- 
tion, i.e., the growth of each column is random and inde- 
pendent of neighboring columns. As a result, the tallest 
column may be next to the shortest column, and w does 
not systematically depend on probed lengthscale. KPZ 
processes can be modeled following the same rules that 
govern Poisson processes, but utilizing "sticky" blocks 
which attach to the first particle they touch. Particles 
falling in short columns can then stick to the side of a 
tall column. This permits short columns that neighbor 
tall columns to grow very quickly; thus, w grows slower in 
a KPZ process than in a Poisson process. Finally, KPZQ 
processes can be modeled following the same rules that 
govern KPZ processes, but with spatial modifications to 
the growth rate. Growth rates are assigned without spa- 
tial correlation such that some regions grow very quickly 
while other regions grow very slowly. This disparity in 
growth rates causes w to increase very quickly. 

To summarize, the growth process of particles de- 
posited at the edge of evaporating suspensions is highly 
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FIG. 4: a. Skewness, 71, multiplied by width, w, plotted 
versus h for e = 1.0, 1.2, 3.5 (solid line, dashed line, dotted 
line, respectively). The line with dashes and dots is jiW = 1, 
the prediction for a Poisson process, i.e., the skewness grows 
linearly with the width. Error bars represent standard devia- 
tion, b. Width, w, plotted versus probed lengthscale, L, for a 
drop containing ellipsoids a = 2.5 with height h = 17 Jim. c. 
Example of the colloidal "Matthew Effect." The trajectory 
of an ellipsoid (e = 2.5) during deposition is shown (red line). 
Initially, it is pushed towards a region that contains very few 
ellipsoids. However, when it adsorbs on the air-water inter- 
face, it is attracted to a particle-rich region, and deposited 
there, d. Dynamic scaling exponent /3 plotted versus particle 
aspect ratio e. Three distinct regimes are readily identified 
(indicated by shading). Predictions for random deposition 
(ft = 0.5), KPZ processes (/? — 1/3), and KPZQ processes 
(/3 — 0.68) are indicated by solid horizontal lines. 



dependent on particle shape. Slightly anisotropic parti- 
cles appear to be deposited by a KPZ process. Spheres 
are deposited by a Poisson like process, until the deposit 
becomes multilayered; at this point the process is best 
characterized as random deposition with surface diffu- 
sion. Finally, very anisotropic ellipsoids induce strong 
capillary attraction on the air-water interface and ap- 
pear to belong to the KPZQ class, i.e., KPZ class with 
quenched disorder. Thus, evaporating drops of suspen- 
sions provide an appealing model system for studying 
theoretical predictions for different growth processes. In- 
teresting future experiments could probe the aspect ra- 
tio boundary regions more comprehensively. Further, the 
various growth processes hold potentially important con- 
sequences for applications involving interfacial colloidal 
aggregates ( e.g., Pickering emulsions and food pro- 
cessing [35j, |36|), and it should be interesting to investi- 
gate how growth processes affect the stability or rigidity 
of colloidal deposits. 
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